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In the light of the recent discovery of a neutron star with a mass accurately determined to be 
almost two solar masses, we carefully examine the maximum mass of a neutron star containing 
hyperons within a relativistic Hartree-Fock treatment. Using the quark-meson coupling model, 
which naturally incorporates hyperons without additional parameters, we find a maximum mass a 
little over 2.1 M©. 
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The recent observation of a 1.97 ± 0.04 Mq millisecond 
pulsar, PSR J1614-2230, by Demorest et al. has set 
the most stringent limit on models of neutron star cores 
so far. This discovery has spurred a re-examination of the 
possibility of exotica such as hyperons, Bose condensates, 
and quark matter playing an important role in models of 
neutron star interiors, owing to a presumed softening of 
the equation of state (EoS) expected in the presence of 
additional degrees of freedom. Historically, this has led 
to expectations of reduced maximum neutron star masses 
for compact objects in hydrostatic equilibrium. 

In this Letter we build on the earlier work of Stone 
et al. Q, who already predicted the existence of neu- 
tron stars containing hyperons with masses as large as 2 
Mq in 2007, to establish a conservative upper limit on 
the maximum mass of such a star. We work within the 
quark-meson coupling (QMC) model [3-5], which has the 
advantage of being derived from the quark level, with a 
very small number of adjustable parameters, while being 
consistent with a broad range of constraints derived from 
hypernuclei as well as normal nuclear properties. We find 
that the stability under variation of the very small num- 
ber of adjustable parameters is such that if a star were 
discovered with a mass significantly above 2.1 Mq, we 
would need to consider more exotic physics, because it 
could not be accommodated within the QMC model. 

We recall that QMC is based upon the self-consistent 
modification of the structure of a baryon embedded in 
nuclear matter. At Hartree level it involves only three 
adjustable parameters which describe the effective cou- 
plings of the a, u) and p mesons to the u and d quarks. 
These are fixed by adjusting them to fit the properties of 
symmetric nuclear matter, namely its saturation density 
and binding energy as well as its symmetry energy. We 
note that the a meson used here simply serves as a con- 
venient representation of the scalar-isoscalar attraction 
arising from two-pion exchange. 

In the most recent development of the QMC model Q , 
the self-consistent inclusion of the gluonic hyperfme inter- 
action led to a very successful description of the binding 



energies of A-hypernuclei — as well as the observed ab- 
sence of medium and heavy mass S-hypernuclei — with 
no additional parameters. We stress that this is achieved 
without any coupling of the strange quark to the a, u) and 
p mesons (which would be OZI suppressed) and with- 
out the need to introduce any further mesons. While 
the model could be supplemented with much heavier 
mesons containing strange quarks Q , Occam's razor sug- 
gests that one should not introduce them if they are not 
needed. 

A clear connection has been established between the 
self-consistent treatment of in-medium hadron structure 
and the existence of many-body [7] or density depen- 
dent [1] effective forces. Dutra et al. Q critically ex- 
amined a variety of phenomenological Skyrme models of 
the effective density dependent nuclear force against the 
most up-to-date empirical constraints. Amongst the few 
percent of the Skyrme forces studied which satisfied all of 
these constraints, the Skyrme model SQMC700, derived 
from the QMC model, was unique in that it incorporated 
the effects of the internal structure of the nucleon and its 
modification in-medium. 

While the earlier study of Stone et al. Q demonstrated 
the importance of exchange (Fock) terms in calculations 
of the EoS of dense baryonic matter in /3-equilibrium, it 
included only the Dirac vector term in the vector-meson- 
nucleon vertices. In this work we include the full vertex 
structure which one might expect to enhance the pressure 
at high density. This is especially so in the case of the p 
meson for which the tensor coupling is much larger than 
that of the u>. Our calculation extends the work of Krein 
et al. [10], who considered nucleons only, by evaluating 
the full exchange terms for all octet baryons and adding 
them in the same way as Stone et al. [2J; as additional 
contributions to the energy density. It also extends and 
complements the important work of Miyatsu et al. who 
calculated the properties of neutron stars within an ear- 
lier version of the QMC model [14]]. In particular, we 
employ the latest version of the model which reproduces 
key hypernuclear properties without the adjustment of 
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coupling constants needed there. Furthermore, we care- 
fully explore the limit on the maximum mass of a neutron 
star containing hyperons while ensuring consistency with 
critical nuclear properties, such as the incompressibility 
of nuclear matter. 

Next we describe the main features of the QMC model 
and explore its parameters arranged into groups ("sce- 
narios") to test the robustness of its predictions. Within 
the QMC model, the baryon energy density, e_B , is given 
by 
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where the effective, in-medium baryon masses, M B , are 
calculated self-consistently for an MIT bag immersed in 
(and in Ref. 0], parameterized as functions of) a mean 
scalar field, designated here by a barred symbol. At a 
given density, a is self-consistently expressed as 
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An additional contribution, 5a, to the scalar field arises if 
we include the Fock terms in the minimization of the en- 
ergy density. This provides only a small correction to the 
mean field, and its effect is included only in the scenario 
denoted 8a. The total hadronic energy density, e B , is the 
sum of baryonic, e B , and mesonic, e CTWp7r contributions, 
for which 
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arise from symmetry considerations, are given in Ref. [2j. 
Note that the ir meson only contributes to the second 
term in Eq. §3§ as it is coupled via a pseudo-vector cur- 
rent. For tuujp, the integrand has the form 
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where A a (k) is the Yukawa propagator for meson a with 
momentum k = p — p'. For the vector mesons, the full 
vertex structure is included in the manner of Ref. llGl as 



2M* 2 v 



(5) 



As usual, the effect of short distance repulsion on the 
Fock terms is simulated by the replacement 
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from which the unit term is subtracted, thus eliminating 
a 5-function. The form factors F" 2 all have the same 
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FIG. 1. Mass vs radius relation for a variety of scenarios de- 
scribed in the text. A BPS [l7j crust has been added in all sce- 
narios. We note that the maximum stable neutron star masses 
(indicated by filled circles) for essentially all full Hartree-Fock 
scenarios are compatible with the observations of Ref. [J (de- 
noted by a black horizontal line with grey dashed error-band. 

dipole form with the cutoff mass A varied from 0.9 to 
1.3 GeV to test the sensitivity. As a further test of the 
model dependence, we consider two choices for the ratios 
of tensor to vector coupling constants n aB = f aB / g aB 
(a G {uj,p}). In scenario kj (consistent with values de- 
rived within QMC) we take these ratios from vector me- 
son dominance (k p n = fpN/gpN = 3.70). Alternatively 
(scenario kii) we take these ratios from the Nijmegen 
potentials (Table VII of Ref. [llfl). with a larger value of 
k p n (= 5.7). 

Of the baryon- meson coupling constants g aB (a), g^B, 
and g pB , only g^ B is density dependent. Its model pa- 
rameterisation [5| is dependent on the free nucleon radius, 
which is taken to be R^ cc = 1.0 fm - with an alternate 
scenario having R l ^ cc = 0.8 fm. The density dependence 
is given by 
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Values of the coupling constants g a N for various 
mesons a and a selection of scenarios considered in this 
work are presented in Table |TJ The couplings g uB and 
g pB are expressed in terms of the quark level couplings 

(8) 



guiB 



g P B = g P N 



9l 



whcr 



represents the number of light quarks in 



baryon B. The a, uj and p couplings to the quarks are 
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Scenario 

hi 
Ku 

A = 1.0 
A = 1.1 
A = 1.2 
A = 1.3 
R = 0.8 
Fock So 
Dirac Only 
Hartree Only 
Nucleon Only 



K 



R 



9aN 9uN 9p (MeV) (km) (M ) (n ) 



10.42 
10.55 
10.74 
11.10 
11.49 
11.93 
11.20 
10.91 
10.12 
10.25 
10.42 



11.02 
11.09 
11.66 
12.33 
13.06 
13.85 
12.01 
11.58 
9.25 
7.95 
11.02 



4.55 
3.36 
4.68 
4.84 
5.03 
5.24 
4.52 
4.52 
7.83 
8.40 
4.55 



298 
299 
305 
312 
319 
329 
300 
285 
294 
283 
298 



12.27 
12.19 
12.45 
12.64 
12.83 
13.02 
12.41 
12.29 
12.56 
11.85 
11.64 



1.93 
1.93 
2.00 
2.07 
2.14 
2.23 
1.98 
1.98 
1.79 
1.54 
2.26 



5.52 
5.62 
5.32 
5.12 
4.92 
4.74 
5.38 
5.5 
5.2 
6.0 
5.82 



TABLE I. Meson-nucleon coupling constants determined for 
our baseline scenario l Ki ' (for which A = 0.9 GeV, and R% ee = 
1.0 fm) and subsequent scenarios in which differences from 
hi are given in column 1. Also shown are the saturation 
incompressibility, K; stellar radius; maximum stellar mass 
and corresponding central density (units no = 0.16 fm -3 ). 

constrained to reproduce a saturation energy per baryon 
of £ sa t = —15.86 MeV and an asymmetry energy coef- 
ficient of a asym = 32.5 MeV at the saturation density 
no = 0.16 fm -3 . The w, p and ir masses are constrained 
to their experimental values, whereas the a mass is taken 
to be 700 MeV. 

For a compact object in /3-equilibrium we solve the 
familiar system of equations for the number densities 
of the baryons and leptons 
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density and pressure are given by the usual formulas 
for a degenerate Fermi gas. In order to obtain the 
neutron star properties shown in Table U we solve the 
Tollmann-Oppenheinmer-Volkov equations for the grav- 
itational mass and radius [IH . The resulting dependence 
of the neutron star mass on radius, for a selection of the 
variations of the model, is shown in Fig. [TJ 



In Table Q] we present the coupling constants, incom- 
pressibilities of symmetric nuclear matter, and stellar 
properties, for a number of variations of the QMC model, 
in each scenario including the a, n, uj and p Fock terms. 
We note that in all scenarios S - hyperons are present in 
significant quantities in the the maximum mass stars. 

It is remarkable that in all of the scenarios investigated, 
the stellar properties are largely consistent, and similar 
to those reported by Stone et al. [2j. Scenarios in which 
the maximum stellar mass lies outside of the range 1.9— 
2.14 Mq correspond to nuclear matter compressibilities 
above the upper limit set in the recent comprehensive 
analysis of giant monopole resonance data [13j. While 
this cannot be true in general, it is certainly the case for 
the QMC model. 

Turning to the effects of the inclusion of the full ex- 
change terms on stellar properties, we find that the 
threshold density for S~ is lowered, while those of A and 
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FIG. 2. Equation of state for a variety of scenarios. Kinks 
occur at significant hyperon threshold densities. The diver- 
gence between the baseline scenario («j) and the 'Hartree 
Only' and 'Dirac Only' scenarios highlights the importance of 
the pN tensor coupling in Hartree-Fock at high density. 




FIG. 3 
(lower) 
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Species fractions (upper) and Fock energy densities 
for our baseline ('ft/') scenario. 



H are raised, as demonstrated in Fig. [3l In all scenar- 
ios there is a greater splitting between the thresholds of 
the S baryons than that found by Stone et al. 0. In 
our baseline scenario ('/tz')> the S~ threshold occurs at 
0.42 fm" 3 , followed by S° at 0.91 fm- 3 . We find that A 
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production is not energetically favoured at densities con- 
sidered here, in agreement with Ref. 14[. Using the Ni- 



jmegen values of tensor coupling strength ('kjj'), the pN 
vector coupling is reduced as the tensor part of the inter- 
action contributes significantly to the symmetry energy. 
The EoS is otherwise largely insensitive to this choice. 
Similarly, it is insensitive to the choice of free nucleon 
radius, despite a moderate impact on the couplings. 

The correction (8a) to the scalar mean field arising 
from the Fock terms decreases the incompressibility by 
13 MeV, yet other observables remain largely unaltered 
by this addition. The cutoff, A, used in form factors 
(which controls the strength of the Fock terms) exhibits 
a more pronounced relationship with the observables in 
Table HI Increasing A beyond 0.9 GeV raises the incom- 
pressibility with the case denoted A = 1.3 GeV already 
exceeding the limit of K < 315 MeV. We stress that 
A could not take a lower value without impacting the 
masses of the a, u, and p. The tt, however, has a much 
lower mass and as such could involve a lower cutoff. We 
investigated this possibility but found only a minimal 
effect on the EoS, as expected from the small contribu- 
tion the pion makes to the EoS at high density (refer to 
Fig. [3]). Overall, increases in the cutoff correlate with in- 
creases in both the saturation incompressibility and max- 
imum stellar mass. 

We stress that the QMC model does not predict the 
appearance of E hyperons at any density where the model 
can be considered realistic. This is in contrast to a num- 
ber of other relativistic models which do predict the E 
threshold to occur, even prior to that of the A 3, 16j |. 
We note that Schaffncr-Biclich [15] considered a phe- 
nomenological modification of the E potential with addi- 
tional repulsion, which significantly raised its threshold 
density. In the case of the QMC model the physical ex- 
planation of the absence of E-hyperons is very natural, 
with the mean scalar field enhancing the repulsive hy- 
perfine force in the bound E. Recall that the hyperfine 
splitting is due to one-gluon-exchange, which determines 
the free E-A mass splitting in the MIT bag model. 

For comparison purposes, we also include a 'Nucleon 
Only' scenario, in which hyperons are artificially ex- 
cluded. In this case the EoS is increasingly stiffer at den- 
sities above 0.4 fin -3 , leading to a large maximum stellar 
mass of 2.26 M Q , consistent with many other nucleon- 
only models. 

It is worth remarking that upon inclusion of the ten- 
sor coupling, the proton fraction increases more rapidly 
as a function of total baryon density. This is likely to 
increase the probability of the direct URCA cooling pro- 
cess in proto-neutron stars. As a further consequence, 
the maximum electron chemical potential is increased in 
this case, which may well influence the production of ir~ 



and K condensates. Changes to the A threshold (occurs 
at higher density with lower maximum species fraction) 
reduces the possibility of H-dibaryon production as con- 
strained by /3-equilibrium of chemical potentials. 

In summary, taking into account the full tensor struc- 
ture of the vector-meson-baryon couplings in a Hartree- 
Fock treatment of the QMC model gives increased pres- 
sure at high density — largely because of the pN tensor 
coupling — while maintaining reasonable values of incom- 
pressibility at saturation density. The conceptual separa- 
tion between the incompressibility at saturation density 
and the amount of pressure or 'stiffness' at higher den- 
sities is critical. It is the latter that leads to neutron 
stars with maximum masses ranging from 1.90 M to 
2.14 Mq, even when allowance is made for the appear- 
ance of hyperons. This suggests that hyperons are very 
likely to play a vital role as consituents of neutron stars. 
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